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The experimental status and theoretical uncertainties of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix describing the charge-changing weak transitions between quarks with charges 
— 1/3 (d, s, b) and 2/3 (u, c, t) are reviewed. Some recent methods of obtaining phases of 
CKM elements are described. 



1 Introduction 

Information about the Cabibbo-Kobayashi-Maskawa (CKM) matrix describing the charge- 
changing weak transitions between quarks with charges —1/3 (d, s, b) and 2/3 (u, c, t) has 
been steadily improving over the years. Despite a wealth of over constraining experiments, no 
significant inconsistencies in its parameters have emerged so far. One seeks greater accuracy in 
the determination of CKM elements not only to expose such inconsistencies, which could signal 
new physics, but also to provide input for an eventual theory of these elements. 
The matrix may be defined in one parametrization^ as 
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0.8, fj ~ 0.36, p ~ 0.19. Many 
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detailed reviews exist! ' ' \ we concentrate on procedures and open questions. 

We shall be concerned with information regarding both magnitudes and phases of CKM 
elements. These are encoded in the angles of the unitarity triangle, illustrated in Fig.^ Current 
fits (not including some CP asymmetries providing information on a and 7) imply la limits^ 
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Thus, although (3 is well known [with the ICHEP 2004 average now (23.3li|)°P, a and 7 are 
more uncertain, ranging over about 40° at the 95% confidence level (c.l.). They can be pinned 
down more precisely using B-B mixing, kaon decays, and CP asymmetries in B decays. 



2 V u d from nuclear, neutron, pion (3 decays 

Our discussion is based on a recent overview^. Nine measurements of nuclear + — > + tran- 
sitions yield an average of \V u d\ = 0.9740(1) (3) (4), where the errors correspond to experiment, 
nuclear theory, and radiative corrections, respectively. Neutron decay gives |V^rf| 2 (l + 3<^)r n = 
(4908 ± 4) s, so using the measured lifetime r n = 885.7(7) s and gA = 1.2720(18) one finds 
\V u d\ = 0.9729(4) (11) (4), where the errors are associated with r n , gA-, and radiative corrections. 
(A very new value ^ r n = 878.5(7)(3) s implies \V ud \ = 0.9757(4) (11) (4).) Pion beta decay 
(tt + — > ?T e + f e ) yields \V u d\ = 0.9739(39); an ongoing experiment at PSI™ seeks to reduce the 
errors further. The overall average (before including the result of Ref.^j) is IV^I = 0.9740(5). 



3 V us : Hyperon and decays; lattice 

Semileptonic hyperon decays, including new measurements of A — * pe~i>, Y,~ — > ne~v, H~ — > 
Ke'u, and H° -> £+e~P, have been analyzed^, with the result \V US \ = 0.2250 ± 0.0027. Sat- 
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Figure 1: Definition of sides and angles of the unitarity triangle. 

isfactory fits to all data have been found without the need for SU(3) breaking, though another 
analysis^^ requires it, obtaining \ V US \ = 0.2199 ± 0.0026. One remaining question is the magni- 
tude of the axial weak- magnetism parameter gi , which is not well constrained by data. 

Several experiments have remeasured K(z decays. Brookhaven E865^^ finds \V US \ = 0.2272=t 
0.0022 rato ±0.0007f.f. ±0.0018j + ( ); radiative corrections were an important part of the analysis. 

Fermilab E832 (the KTeV Collaboration) 13 obtains \V US \ = 0.2252 ± 0.0008 K TeV ± 0.0021 ext , 
where "ext" refers to all errors external t o th ose in KTeV, such as uncertainties in form factors, 
making use of the radiative corrections in^l 

The value of \V ud \ = 0.9740 ± 0.0005, when combined with the Particle Data Group El 
value \V US \ = 0.2200 ± 0.0026 agreed poorly with CKM unitarity: \V ud \ 2 + \V US \ 2 + \V ub \ 2 = 
0.9971±0.0015. while the expected value of | V us | from unitarity is | V us \ = (1— |Kd| 2 — iKfel 2 ) 1 ^ 2 = 
0.2265±0.0023. The above two JQ3 results, as well as thatP^from hyperon decays (see, however, 
Ref.El), are much more consistent with unitarity. So, too, are very recently reported results 
from the KLOE detector at the DA$NE e + e~ collider at Frascati 15 . The NA48 Collaboration 
^reports a value \V US \ = 0.2187 ± 0.0028 using a value of /+(0) about 1.35% higher than that 
used by KTeV, so the situation is still not entirely settled. 



4 V cd and V cs : Charm and W decays; neutrino production 

The CLEO III detector has reported rates and spectra for D° — » ir~l + v and D° — > K~£ + iM^ 
implying \f« (0)| 2 | V cd | 2 /|/_f(0)| 2 | V cs \ 2 = 0.038±g;gg?±g;ggi . Progress in n f = 3 lattice QCD by 
the Fermilab-MILC CollaborationQSlhas yielded form factors f£~™(0) = 0.64(3)(5), f^ K (0) = 
0.73(3) (6), where the lattice errors are statistical and systematic, leading when combined with 
previous measurements of nonstrange and strange charm decays to \V C( i\ = 0.239(10)(19)(20), 
\V CS \ = 0.969(39) (78) (24), where the final errors are experimental. 

Charm production by neutrinos (signaled by dileptons) leads^l to £>(c — > i + X)|V C( i| 2 = 
(4.63 ± 0.34) x 10~ 3 so that with B(c -> i + X) = (9.23 ± 0.73)% one has \V cd \ = 0.224 ± 
0.012. New measurements by the CLEO-c Collaboration of B(D° -> [K~,TT-]e + v) and B{D+ -> 
[K*°, p a ]e + v) with 57 pb _1 at the ^"(3770) rep rese nt a further source of information on 
V c d and V C s when combined with the lattice results ^1 For example, the CLEO-c result for 
B(D° -> irev)/B(D° Kev) is 0.070 ± 0.007 ± 0.003, to be compared with CLEO-iU's 0.082 ± 
0.006 ±0.005. Higher-luminosity running in e + e~ — > ip"(3770) — > DD is to begin in September, 
and further improvements are envisioned. The CLEO-c Collaboration eventually hopes for 3 
fb _1 at the ifi" . Lattice errors will be the limiting factor in extracting \V cd \ to 0(1%). 

Charm-tagged W decays at LEPII (ALEPH, DELPHI) - have given the value \V CS \ = 0.97 ± 
0.09 ± 0.07. Measurement of the leptonic branching ratio B(W — ► iv) and the assumption of a 



standard pattern of W decays can be used to improve this estimate through the relation 



^Z\ v ij?j {i = u,c; j = d,s,b) (3) 

which implies l^ijl 2 = 2.039 ± 0.025 and hence \V CS \ = 0.996 ± 0.013 when contributions of 
other CKM elements are subtracted. The study of Cabibbo-favored (c — * s) decays of charmed 
particles at CLEO-c will provide \V CS \ at 0(1%) accuracy if the C goal is achieved and if lattice 
gauge theory continues to progress. 

At this juncture one can't see a violation of two-family unitarity: 

\Vud? + \Vus? = \V cd \ 2 + \V CS \ 2 = \V ud \ 2 + \V cd \ 2 = \V US \ 2 + \V CS \ 2 = 1 . (4) 

This reflects on the very hierarchical structure of the CKM matrix. It implies that violations of 
unitarity may be too small to signal the presence of additional quark families. 
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5 V c b' b — > c inclusive and exclusive decays 

At the quark level, the b — > c semileptonic decay rate is simple (neglecting m^): 

r(6^ c ^-) = I ^|F cfe | 2 /f^J , /(x)^l-8x + 8x 3 -x 4 -12x 2 lnx . (5) 

However, initial and final states contain real hadrons, so one must employ perturbative QCD; 
expansions in inverse heavy quark mas s ( " HQE" ) ; and moments in lepton energy, hadron mass, 
and photon energy in b — > yy (se e, e.g. ,1=^ for a review of these techniques) to infer the inclusive 
value \V cb \ = 0.0421 ± 0.0013^ which is quoted in 5 as 0.0420 ± 0.0006 sta t ± 0.0008 theo . More 
recent contributions include \V c b\ = (41.4 ± 0.4 sta t + 0.4hqe ±0.6 thco ) x HT^and \V cb \ = 
(42.4 ± 0.8 stat+H QE ± (> 0.8) theo ) x W~ 3 M 

The best source of \ V c f,\ from exclusive decays is the process B — > D*lv£, whose rate is 

2 



fa = -^\Vcb\ 2 (m B - m* D ) 2 m 3 D *V^ig(w)\T(w)\ 2 . (6) 

One Isgur- Wise form factor ^(w) ~ J r (\)[l+p 2 {w — 1) . . .], a function of the variable w = vb-vd* , 
governs the process; phase space is described by a function Q(w) with G(l) = 1. The form factor 
at to = 1 is given by = t)qcd\X + 0(l/ml)] = 0. 913^035 as calculated using lattice and 
HQET estimates^. 

A compilation of values of .F(l)| V^l^ is shown in Fig. |5J The latest average is .F(1)|T4&| 2 = 
(37.7 ± 0.9) x 10~ 3 , with the form factor slope parameter p 2 = 1.56 ± 0.14, leading to \V cb \ = 
(41.4 ± l.Oexp + l-8th) x 10 -3 . This form factor shape is consistent with the rates for B° — > 
/)(*)- 7r + ) I)(*) _ Z)^*- >+ estimated with a factorization approach It is consistent with the 
inclusive \ V c b\ value. Although it currently has larger errors, progress in experiment and lattice 
estimates may eventually make this the best source of \V cb \. 



6 V u b- B — » u inclusive and exclusive decays 

The semileptonic decay b — * uvgl~ is the source of information on V ub from inclusive decays. 
However, T(b — > uD^i~) is only about 2% of T(b — > cug£~). Several strategies have been used to 
extract the b — ► uvt>l~ contribution^^, including measurement of leptons beyond the b — » cvil~ 
end point, reconstruction of hadronic masses Mx < Mp, cutting on q 2 = m 2 u , and, most 




Figure 2: Values of JT(1)|V C (,| vs. form factor parameter p 2 measured by various groups. 

recently, light-cone-inspired kinematic cuts^A As in b —* ci>i£~, the b — * S7 photon spectrum 
helps to pin down hadronic uncertainties. 

A compilation of values of | V^,|^is shown in Fig. |31 The bottom 5 points lead to an average 
\V u b\ = (4.66 ± 0.43) x 1CP 3 . This is notably higher than the value from exclusive b — > u decays 
(e.g., B — > iruii). The top 6 plotted points give an average of \V u b\ = (3.26 ± 0.62) x 10~ 3 from 
exclusive decays. In this average I assigned an overall systematic error of ±0.60. Combining 
with the inclusive value and including a scale factor^! S = 1.86 = v5?> I find \V u b\ = (4.21 ± 
0.66) x 10~ 3 . The slight discrepancy between inclusive and exclusive values merits caution; 
quark-hadron duality may not be valid if extended down to hadronic masses which are so low 
as to be represented by discrete states like tt and p. 

The measurement of the spectrum for B° — > 7r~Z + z^, e.g., as performed by CLEO^, can 
be useful not only in extracting \V u b\ from lattice gauge theory for m fa ctor results (typically 
obtained for high q 2 ), but also in testing factorization in B° — > 7r + 7r _ l2ni. 

7 V t d- B°—B° mixing; progress on decay constants 

Loop diagrams with quarks i, j = u,c,t in the intermediate state allow bd <-> db transitions at 2nd 
order in weak interactions. The t quark dominates, so this mixing provides information on \ Vtd\- 
The predicted splitting Am^ between mass eigenstates in the B°-B system is proportional to 
a parameter fg (the B meson decay constant) governing the matrix element of the bd <-> db 
operator between physical meson states, and to a parameter Bb equal to 1 if W exchange 
diagrams dominate. A lattice estimate fsV^B = (228 ± 30 ± 10) MeV and the experimental 
value Arrid ~ 0.5 ps _1 imply the 95% c.l. range - \V t d\ = (8.26jz^'yg) x 10~ 3 , equivalent to 
|1— p— ifj\ = 0.891q2o- (The upper limit on these quantities is governed by a lower limit on B s -B s 
mixing to be discussed below.) For comparison, the 95% c.l. range^of \V u b\ = (3.87lQg 3 ) x 10~ 3 
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Figure 3: Values of \V u b\ measured by various groups. 



implies = 0.40t 06- Information from B(B — > p~/)/B(B — > K*j) also constrains |V^/T4 S |. 

One expects, for example^ B(B° p°j) = (0.64 ± 0.23) x 10 -6 . A recent report^ finds 
B(B° -* p°7) < 0.4 x 10" 6 and \V td /V ts \ < 0.19 at 90% c.l. 

The new CLEO value f D = (201 ± 41 ± 17) MeV^ has an accuracy approaching that of 
lattice calculations. With increased integrated CLEO luminosity, this quantity will be measured 
precisely enough to test those calculations, lending weight to their predictions for fs- 

As mentioned, \Vtd\ is quoted with an asymmetric error since it is restricted on the positive 
side by B s -B s mixing. This quantity is described by the same diagram as B°-B mixing with 
the substitution d — > s. We assume Vt s — —V c b- Then the lower limit Am s > 14.5 ps _1 and the 
SU(3)-breaking estimate £ = fB s \JBsJ JbVBb = 1-21 ± 0.06 implies a lower limit on Vt s /V t d 
and hence an upper limit on Vtd- The (±2cr) prediction of Ref.^is Am s = 17.8l2 5 7 2 P s_1 - 

The mixing of B s and B s can pro ceed via on-shell shared intermediate states, as described 
in Fig. One expects (see, e.g.,'""' and references therein) Ar s ~ — Am s /200 (~ m?/mf) 
or Ar s /f s ~ 0.18(/b s /200 MeV) 2 in lowest order, but a state-of-the-art calculation^ finds 
Ar s /f s = 0.12 ±0.05 for f Bs = 245 MeV. The CDF Collaboration has recently reported a value 
of 0.65±o || ±0.01^1. Since f Bs is expected to be no larger than about 300 MeV, the CDF result 
is somewhat larger than anticipated, but not yet in serious conflict with theory. 



Figure 4: One graph describing B s ~B a mixing, which can correspond to a real intermediate cess state. 
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Figure 5: Summary of 95% c.l. constraints on (p, fj) due to \V u b\ (short-dashed circles), B°-B° and B s -B 3 mixing 
(dashed circles), CP-violating K°-K mixing (solid hyperbolae, ±1ct limits), and sin2/3 from CP asymmetry in 

B° — * J/ipKs and related processes (dash-dotted lines). 

8 Constraints from Kl, K + decays 

CP-violating K°-K° mixing is dominated by top quarks in second-order-weak loop diagrams. 
The parameter ek describing this mixing depends mainly on Im(V^) and hence measures ap- 
proximately fj(l — p), with charmed quarks supplying a small correction. Convenient expres- 
sions in Ref. imply that (p, fj) lies between the la boundaries ^(1.38 — p) = 0.28 and 
7/(1.26 — p) = 0.50. These constraints and those on (p, fj) from \V u b\, \Vtd\, and sin 2(3 select a 
region around (0.19±g;^ ,0. 36^^04), shown as the plotted point in Fig. [5J 

Note the consistency of all these determinations. The main uncertainty is associated with 
the value of |1 — p — irj\. Direct CP violation in neutral kaon decay (governed by a parametere'/e) 
is seen but provides no useful constraint as a result of hadronic uncertainties. 

In K — > ttvu decays, higher-order weak diagrams govern the weak quark transition s — > dvv. 
For K + — > tx^vv the top quark in the loop dominates but there is also a charm contribution, 
so th at th e rate measures 1 1 . 3 — p — ifj\. An experiment (E787, E949) at Brookhaven sees three 
events^^, corresponding to a branching ratio B(K + — > it + vi>) 

= (1-47±J;1§) x 10- 10 . Comparing 
with the Standard Model prediction^ of (0.78 ±0.12) x 10~ 10 , the result is still consistent with 
a large (p, fj) region. Proposals for Fermilab and JPARC seek a sample of 100 events which 
could determine |1.3 — p — ifj\ to 5%. 

The process Kl — > -k^vv is purely CP-violating and measures ff. One expects B = (3.0 ± 
0.6) x 10 _11 l"2. An experiment at KEK (PS E391) has taken data whose single-event sensitivity 




Figure 6: Left: B° — » ty°Ks; solid ellipse: allowed region boundary. Right: B° — > r( Ks\ solid curve: bounds 
based on flavor SU(3) and measurements of other processes; dot-dashed curve: bounds prior to latest set of such 
measurements; dashed curve: neglecting processes involving spectator quark. 



should be an order of magnitude above the Standard Model valueM*", this represents tremendous 
progress in the past few years. The eventual goal of proposed experiments at JPARC and 
Brookhaven is to be sensitive at the Standard Model level. 



9 Vts and Vtb' b — > 57, top quark decays, unitarity 



One may obtain a lower limit on |V^T4 S | from B s -B s mixing, the assumption that £ = 1.2 ±0.1 

* 
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and Bd-Bd mixing, yieldingly | V^^VtsI > 0.034. An upper limit can be extracted from the top 



quark contribution to b — > 57^! V#JVf S = —0.047 ± 0.008, which has the expected sign. 

The fact that top quark decays t — > M + V£ dominate over those with no b has been used by the 
CDF Collaboration to conclude on the basis of Run I data that | Vtb| 2 /(|^4d| 2 + \ Vt s \ 2 + \ Vtb\ 2 ) = 
0.94^Q24- A slightly less stringent limit has been placed for a sample of 108 pb _1 of Run II 
data: 0.54+g;^. The assumption that CKM matrix is 3 x 3 is crucial to interpretation of these 
results as a useful bound on |Vtb| 2 . The lower bounds on magnitudes of the third row of the 
CKM matrix are almost non-existent if there are more than 3 families 4 . 



10 Processes giving j3 

Many experiments measure a time-dependent CP asymmetry in B decays of the form A(t) = 
— C cos(Ami) + Ssin(Ami). Results involving the subprocess b — > ccs, such as B° — > J/ipKs, 
show beautiful agreement with Standard Model fits, as illustrated in Fig. \5\ 

For processes dominated by the b — ► s penguin amplitudes one expects C = 0, S = sin(2/?) = 
0.74 ± 0.05. One can estimate contributions from other amplitudes usi ng new measurements 
and flavor SU(3), as in the processes B° — > Tr ^^'^! anc [ go _ ¥ ^'^gEDJ, We show in Fig. 
El examples of bounds on allowed deviations from the Stand ard Model predictions. We have 
applied a scale factor in averaging rj 'K s results from BaBarEU (S v , Ks = 0.27 ± 0.14 ± 0.03, 
C n , Ks = -0.21±0.10±0.03) and BelleES (S n > Ks = 0.65±0.18±0.04, Cr,'K. = 0.19±0.11±0.05), 
so the plotted data points have larger error bars than quoted in, e.g., Ref.'"'). The corresponding 
values for tt°K s are Sj> Kg = 0.35±g;|g ± 0.04, ^0^ = -0.21 ±0.10 ±0.03 (BaBarED) S n o Ks = 
0.30±0.59±0.11, C n o K = 0.12±0.20±0.07 (Belle 42 ). There is not, in my opinion, evidence yet 
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Figure 7: Results of % 2 fits to B — > PP (left) and B —* VP (right) decays. In PP fit the lower curve corresponds 
to a fit with two extra parameters. In VP fit the solid curve corresponds to a fit with no constraint between two 
independent penguin amplitudes, while the dot-dashed and dashed curves correspond respectively to constraining 
these amplitudes to be relatively real, and equal and opposite. 



for substantial deviations from the Standard Model predictions, but the situation bears watching 
both in the processes illustrated in Fig. Eland in the decay B° — ► fiKg, also dominated by the 
b — > s penguin. With improved data one could well have evidence for new physics. 

11 B° — > 7r + 7r - , ir^pT, p + p~ and a 

The time-dependent parameters S in the processes B° — > 7r + 7r~, B° — » p^ir^, and B° — ► p + p~ 
would measure just sin 2a if one could neglect the effect of penguin "pollution" of the dominant 
tree amplitudes. In B — > -kit the penguin/tree amplitude ratio was estimated to be ~ 70% from 
B — > i^7r decays^} leading to the estimate a = (103 ±17)°. In P° — ► p^vr^ the relative penguin 
contribution is found to be less^l leading to a = (95±16)°, while i n B ° — > p + p _ a recent BaBar 
measurement, combined with an estimate of penguin effects, gives'^ a = (96 ± 10 db 4 db 13)°. 

12 B -» PP, VP and 7 

One can perform fits to rates and CP asymmetries in flavor SU(3) based on amplitudes des- 
ignated by T, P, C, . . . (denoting tree, penguin, color-suppressed, etc.), for B — > PP^Hl and 
P — > ypelJ decays, where P, V denote light (pseudoscalar, vector) mesons composed of u, d, s 
quarks. The B — > PP fits involve 26 observables while the P — > yP fits involve 34. The results 
are shown in Fig. [7| We place greatest reliance on the fits dscribed by the solid curves, corre- 
sponding to Xmin/ d - f - = 16.0/13 (left) and Xmin/ d - f - = 20.5/22 (right). The fit to PP decays 
gives a rather shallow minimum around 7 = 54° whose stability we do not yet trust, while the 
fit to VP decays gives 7 = (63 ± 6)°. The combined fit gives lo cal minima at 7 ~ 55° and 61° 
after being updated in view of new data presented at ICHEP 04^3. 

13 Beyond the 3 x 3 CKM matrix 

Does a fourth family of quarks and leptons exist? Any fourth neutrino must be heavy; only 3 
light neutrinos are seen in Z decayESI A direct search for b' heavier than Z (looking for b' — > bZ) 
by CDF at the TevatronEOl excludes 100 < m(b') < 199 GeV/c 2 . However, looking outside the 
familiar pattern, existing quark-lepton families belong to 16-dimensional multiplets of the grand 



unified group SO(10), consisting of 1 + 5* + 10 of SU(5). The smallest representation 27 of ~Eq, 
an interesting group sontaining SO(10), is 16 + 10 + 1 of SO(10). The 10 of SO(10) consists 
of isosinglet quarks ll h,h v with charge Q = ±1/3, and isodoublet leptons. The SO(10) singlets 
are candidates for sterile neutrinos, one for each family. 

The exotic h quarks can mix with b and push its mass down wit h re spect to Production 
signatures of hh at the Tevatron and LHC have been investigated ^ . Through decays of h to 
Z + b, W + t, and possibly Higgs + b, one can reach h masses 270-320 GeV/c 2 at the Tevatron. 



14 Summary 

(1) The CKM unitarity relation \ V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1 seems valid, though questions remain. 

(2) Improved V c d and V cs values still are consistent with V c d = —V us and V cs = V U( i- 

(3) Errors on V& ~ 42 x 10~ 3 from inclusive analyses are being reduced; a la error of about 3% 
is a conservative guess. Lattice QCD will help exclusive analyses catch up. 

(4) Lower V ub values [(3.26 ± 0.62) x 10~ 3 ] are obtained in exclusive b — > u decays than those 
[(4.66 ± 0.43) x 10~ 3 ] in inclusive decays. One needs a better understanding of form factors and 
quark-hadron duality. An average ~ 4.2 x 10 -3 is known at present to ~ 15%. 

(5) Errors in \V t( i\ — (8.3tJ'g) x 10~ 3 (95% c.l.) are due mainly to lattice errors in /bV^b- One 
expects a major improvement from detection of B s -B s mixing. 

(6) The results of an SU(3) fit giving 7 = (63± 6)° are being validated by the most recent 
meas ure ments of rare B decay branching ratios^l The detection of b — > d penguins at expected 
rates^^l has provided partial confirmation. 

(7) There is no evidence against Vt s ^ — V c b, Vfb — 1. 

(8) It's time for a theory of quark masses and CKM elements! 
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